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ARTICLE INFO 


ABSTRACT 


Article history: 

Received 26 June 2013 
Received in revised form 
28 August 2013 
Accepted 31 August 2013 


Keywords: 

Ozone Monitoring Instrument 
Nitrogen dioxide 
Air quality 

Space-based observations 
Power plant emissions 
Emission control devices 


We show that Aura Ozone Monitoring Instrument (OMI) nitrogen dioxide (N0 2 ) tropospheric column 
data may be used to assess changes of the emissions of nitrogen oxides (NO x ) from power plants in the 
United States, though careful interpretation of the data is necessary. There is a clear response for OMI 
N0 2 data to NO x emission reductions from power plants associated with the implementation of 
mandated emission control devices (ECDs) over the OMI record (2005-2011). This response is scalar for 
all intents and purposes, whether the reduction is rapid or incremental over several years. However, it is 
variable among the power plants, even for those with the greatest absolute decrease in emissions. We 
document the primary causes of this variability, presenting case examples for specific power plants. 

Published by Elsevier Ltd. 


1. Introduction 

In response to federal and state regulations, total emissions of 
nitrogen oxides (NO* = NO + N0 2 ) decreased since the late 1990s 
by 47% in the United States (US). Emissions from electric power 
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generation and highway vehicles, two of the largest sources, 
decreased by 68% and 43%, respectively (http://www.epa.gov/ttn/ 
chief/trends/index.html). The US Environmental Protection 
Agency (EPA) issued the 1998 NO* State Implementation Plan (SIP) 
Call with the intent to reduce emissions in 22 eastern states during 
the summer season so as to decrease ozone. In 2005, it issued the 
Clean Air Interstate Rule (CAIR) for 27 eastern states with the goal 
to decrease NO x emissions even further from power plants. Indi- 
vidual state rules and court orders have also contributed to power 
plant emission reductions. The mobile source of NO x emissions has 
declined nationwide as a result of the requirements of the Clean Air 
Act Amendments (CAAA) of 1990, specifically the Tier 1 (phased-in 
between 1994 and 1997) and more stringent Tier 2 (phased-in 
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between 2004 and 2009) standards, and the gradual turnover of 
the fleet of light-duty vehicles (e.g., Dallmann and Harley, 2010; 
McDonald et al., 2012). 

Satellite observations confirm that NO 2 columns over power 
plants and urban areas in the US have declined as a result. Kim et al. 
(2006) used both the European Remote-sensing Satellite-2 (ERS-2) 
Global Ozone Monitoring Experiment (GOME) and Envisat Scan- 
ning Imaging Absorption spectroMeter for Atmospheric CHartog- 
raphY (SC1AMACHY) N0 2 column data to infer that NO* emissions 
from power plants in the Ohio River Valley decreased from 1997 to 
2005 by about 35%, which is consistent with reported emission 
changes from the Continuous Emissions Monitoring System 
(CEMS). Kim et al. (2009) were the first to show that N0 2 columns 
from a model of chemistry and transport (CTM) using CEMS data 
were consistent with columns from three retrievals (i.e., the Uni- 
versity of Bremen Ozone Monitoring Instrument (OMI), National 
Aeronautics and Space Administration (NASA) OMI operational 
product, and SCIAMACHY) over 13 isolated power plants in the 
western US in 2005. Russell et al. (2012) used N0 2 column data 
from the OMI Berkeley High Resolution (BEHR) retrieval algorithm 
to infer that NO* emissions changes from large power plants were 
variable because of regionally-specific regulations, decreasing by 
26 ± 12% from 2005 to 2011. They estimated an average total 
reduction of 32 ± 7% in N0 2 for US cities from 2005 to 2011 with a 
34% decrease in N0 2 from mobile sources. They attributed part of 
the observed decline to the turnover in the mobile source fleet and 
part to the global economic recession that began in 2008. 

To comply with federal and state requirements, emission control 
devices (ECDs) were installed on power plants, which create a 
natural experiment to assess the response of the satellite-observed 
tropospheric N0 2 column to a known, and oftentimes rapid and 
significant, change in a power plant’s emissions. For instance, in 
Selective Catalytic Reduction systems (SCRs), ammonia is mixed 
with the flue gas before entering the reactor so that ammonia and 
NO* react to form nitrogen and water. Other techniques to reduce 
NO* emissions include the installation of Low NO* Burners (LNBs) 
and Rotating Opposed Fire Air (ROFA) devices, which may be used 
in combination with SCRs. ECDs remove up to 90% of NO* from the 
effluent. 

The purpose of this study is to use Aura OMI data (2005—2011 ) 
to understand the response of the N0 2 column to a change in a 
power plant’s emissions; hereafter, we refer to this as the 
“Response”. As we will show, the Response is scalar, as the change 
in the column is a linear function of the change in emissions for all 
intents and purposes. However, there are variations in the magni- 
tudes of the Responses. We document the primary sources of these 
variations. Quantifying the Response and understanding the pri- 
mary drivers of its variability for power plants in the US will allow 
for 1 ) confidence in the assessment of the impact of ECDs on air 
quality and 2) better estimation of NO* emissions from large point 
sources in other regions of the world where estimates of emissions 
are often highly uncertain. 

2. Data and method 

2.1. OMI NO 2 column data 

The OMI is on board the Aura satellite, which was launched on 
July 15, 2004 into a sun-synchronous polar orbit. It measures direct 
and backscattered solar radiation in the UV— visible range from 264 
to 504 nm (Levelt et al., 2006) and provides early afternoon (local 
time 1300—1430) N0 2 columns at a spatial resolution of up to 
13 x 24 km 2 with global coverage within two days. We use the OMI 
operational tropospheric N0 2 column data product (version 2.1, 
collection 3) from 2005 to 2011, which is available from the NASA 


Goddard Earth Sciences, Data and Information Services Center (GES 
DISC; http://disc.sci.gsfc.nasa.gov). The early releases of the two 
main OMI products of N0 2 , one from NASA and the other from the 
Royal Netherlands Meteorological Institute (KNMI), showed large 
differences for some regions (Lamsal et al., 2010). This current 
version represents substantial OMI retrieval algorithm improve- 
ments (Boersma et al., 2011; Bucsela et al., 2013, and references 
therein) from its preceding version 1 . 0 , so that it is now feasible to 
derive quantitative information about NO* emissions from large 
point sources (Streets et al., 2013). The current, refined retrieval 
algorithms of both research groups, though different in their ap- 
proaches, now produce very similar columns. 

Retrieval of tropospheric N0 2 columns involves (1) retrieval of 
N0 2 abundance along the viewing path (slant column) with a Dif- 
ferential Optical Absorption Spectroscopy (DOAS) fit (Platt, 1994) in 
the 405—465 nm wavelength range, (2) computation of an air mass 
factor (AMF) by integrating the relative vertical distribution (shape 
factors) of N0 2 weighted by altitude-dependent scattering weights 
for N0 2 (Palmer et al., 2001), (3) removal of cross-track artifacts 
(stripes) resulting from insufficient calibration in the OMI back- 
scattered reflectances, and (4) separation of stratospheric and 
tropospheric N0 2 components (Bucsela et al., 2013). 

The tropospheric AMF is sensitive to the a priori N0 2 profile 
shape. The retrieval of the operational N0 2 product uses N0 2 shape 
factors generated from the NASA Global Modeling Initiative (GMI; 
http://gmi.gsfc.nasa.gov/) CTM at 2.5° longitude x 2° latitude res- 
olution grids. In this work, we use the N0 2 product discussed in 
Lamsal et al. (2013) that was generated with high resolution (0.67° 
longitude x 0.5° latitude) over the US. N0 2 shape factors were 
derived from a nested-grid GEOS-Chem CTM (http://acmg.seas. 
harvard.edu/geos/) simulation and scattering weights for N0 2 . 
Use of N0 2 shape factors from the nested simulation improves the 
representation of vertical distributions, including those of the 
elevated plumes of power plants (Lamsal et al., 2013). The errors in 
the individual pixel tropospheric N0 2 columns under clear-sky 
conditions are estimated to be 30% (Boersma et al., 2004). 

The OMI tropospheric N0 2 columns agree with in situ and 
ground-based measurements within 20% (Lamsal et al., 2013; 
Bucsela et al., 2013). Individual clear-sky (i.e., cloud fraction < 0.3) 
data not affected by the so called “row anomaly” (Dobber and Braak, 
2010) were allocated by area-weights into 0.1° longitude x 0.1° 
latitude grids. The row anomaly is the result of a partial blockage of 
the field of view of the OMI, which lengthens the time necessary to 
obtain global coverage from one day to two days. For consistency 
over the OMI record, we restricted our analysis to scan positions 10— 
23 as they are unaffected by the row anomaly. 

There are rarely ideal conditions for assessing emissions 
changes from power plants from space, so we used all available 
data, regardless of season. Lu and Streets (2012), and references 
therein, recommend using data only for summer (e.g., the policy- 
relevant ozone season of May— September) for a variety of rea- 
sons. For instance, the chemical lifetime of NO* tends to be shortest 
in summer, which has the advantage that a facility’s emissions are 
convolved less with NO* from other sources than in other seasons. 
However, there are disadvantages to using data only for the ozone 
season, such as the stratospheric contribution to the total N0 2 
column, and the associated error, is seasonally greatest and 
important. 

For the purposes of this study, it was not practical to restrict our 
analysis to the ozone season as many facilities, including some of 
the largest emitters, were already operating ECDs during the ozone 
season at the start of our study period, especially in the eastern US. 
In Section 3, we show that the relationship between a facility’s NO* 
emissions and the OMI N0 2 column over the power plant is much 
stronger in the southern US than in the northern US, where the 


104 


B.N. Duncan et al. / Atmospheric Environment 81 (2013) 102-111 


seasonal variation in the chemical lifetime of NO* is more pro- 
nounced. We discuss the implications of using all available data in 
Section 3.2.5. 

2.2. Selection of power plants 

We identified the top 100 highest-emitting power plants in 
2005 based on the US national NO* emissions inventory (http:// 
www.epa.gov/ttn/chief/eiinformation.html). If there is more than 
one power plant within a 0.4° longitude x 0.4° latitude gridbox, we 
combined and treated them as one facility. Then, we used the 
Emission Database for Global Atmospheric Research version 4 
(EDGAR v4; http://edgar.jrc.ec.europa.eu/), which is for 2005 and is 
available on a 0.1° longitude x 0.1° latitude resolution grid, to select 
the power plants least affected by other industrial sources within a 
0.4° longitude x 0.4° latitude area around the facility. Lu and 
Streets (2012) found that the agreement between N0 2 columns 
and NO* emissions improves with increasing relative contribution 
of the power plant’s emissions to total NO* emissions from all 
sources within a N0 2 column gridbox; they refer to this quantity as 
fpower- For this study, we required that/ powe r > 0.90, which elimi- 
nated 45 of the top 100 highest-emitting power plants from our 
analysis. The locations of the facilities used in this analysis are 
shown in Fig. 1. Information for each power plant is given in Table 1. 

The characteristics of a plume from a power plant depend on 
variations in meteorology (e.g., “plume meandering”; Beirle et al., 
2011) and chemistiy, so we used the maximum value of the 
plume for each overpass whether in the gridbox containing the 
facility or in adjacent gridboxes (i.e„ a 0.3° longitude x 0.3° latitude 
area) to create monthly averages. We found that the correlation 
between the change in emissions from a power plant and the 
concomitant change in the OMI N0 2 column is best for this fine grid 
resolution (i.e., 0.1° longitude x 0.1° latitude) as compared to more 
coarse resolutions (e.g., 0.25° longitude x 0.25° latitude). 

2.3. Definition of the response (p) 

In order to reflect the relative contributions of NO* emissions 
from a power plant to the total NO* emissions (and column), we 
define the following parameters: E T , E pp , and E°, which represent 


NO* emissions from all sources within a gridbox, the power plant, 
and sources other than the power plant, respectively, where 
E r = E pp + E°. Similarly, N0 2 , NO pp and N0 2 represent, respectively, 
the total N0 2 column within a gridbox, the portion of the column 
associated with the power plant, and the portion of the column 
associated with all other sources, including N0 2 advected into the 
gridbox, where N0 2 = N0 2 p + N0 2 . 

As defined in the introduction, the Response (p) is given by: 

p = ANO£ p /AE pp = (ANOj - AN0 2 ) j (AE r - AE°) (1) 

where A represents the change in N0 2 column or NO* emissions. 
Rearranging Equation (1) into linear form: 

ANOj = p*(AE r - AE°) + AN0 2 (2) 

where p is the slope of the line and AN0 2 is the y-intercept. In the 
ideal situation where AE° and AN0 2 equal zero, Equation (2) sim- 
plifies to: 

ANO£ = p*AE pp (3) 

In this case, AN0 2 , such as determined from OMI data, is solely 
due to AE PP . AE PP in all figures and Table 1 in this manuscript is the 
sum of emissions for those days where OMI data are available. 
Consequently, AE PP is less than the total change in a facility’s 
emissions over a given time period. 

Though Equations (1)— (3) are rather straightforward, p is a 
complicated parameter that is a function of the chemical lifetime of 
NO*, meteorology, and the factors that affect the partitioning of NO* 
into NO and N0 2 (e.g., Martin et al., 2003; Stavrakou et al., 2008; 
Beirle et al., 2011; Lamsal et al., 2011; Walter et al., 2012; Zhou 
et al., 2012). The dependence of p on these factors is discussed in 
the Supplemental Material. Accounting for the complexities of p re- 
quires a CTM, ideally with a plume-in-grid technique, to properly 
treat the evolution of a power plant's plume. In Section 3, we show 
that this onerous step is not necessary for our practical application, 
particularly given the large uncertainties associated with the OMI 
data discussed in Sections 2.1 and 3.2. In practice, p is relatively stable 



Fig. 1 . The locations of the facilities listed in Table 1 with the facility identifier beside each point. The magnitude of AE PP (kTon) is indicated by the size of the circle around the 
square indicating the facility location. The color of each square corresponds to the correlation (r 2 ) shown in Table 1. The horizontal dashed line indicates 36.5°N latitude, which is the 
boundary between the southern and northern US in Fig. 3. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


B.N. Duncan et al. / Atmospheric Environment 81 (2013) 102-111 


105 


Table 1 

Information for individual power plants. 


ID 

Facility name 

State 

Lati- 

tude 

Long- 

itude 

AE PP 

(kTon) 

AE PP 

m 

ANO 2 (xlO 15 
molec cm -2 ) 

ANO 2 

(%) 

p (xlO 15 molec 
(cm 2 kTon) -1 ) 

r 2 

Std 

Err 

N <3 

1 

Crystal River 

FL 

29.0 

-82.7 

5.4 

69.5 

1.5 

36.6 

0.28 

0.91 

0.23 

0 

2 

Cardinal + WH Sammis 

OH 

40.4 

-80.6 

4.9 

69.9 

2.6 

35.5 

0.44 

0.77 

0.70 

15 

3 

Paradise 

KY 

37.3 

-87.0 

4.9 

69.0 

2.0 

36.1 

0.31 

0.71 

0.58 

2 

4 

NewMadrid 

MO 

36.5 

-89.6 

4.5 

74.7 

1.1 

26.5 

0.24 

0.85 

0.26 

2 

5 

Jeffrey EC 

KS 

39.3 

-96.1 

2.9 

44.3 

0.4 

11.4 

0.12 

0.82 

0.23 

2 

6 

Gorgas + James H. Miller Jr. 

AL 

33.6 

-87.1 

3.9 

49.9 

1.3 

28.1 

0.34 

0.91 

0.29 

2 

7 

Monroe 

MI 

41.9 

-83.3 

3.8 

59.7 

2.0 

23.4 

0.48 

0.58 

0.83 

14 

8 

Bowen 

GA 

34.1 

-84.9 

4.1 

69.8 

1.8 

30.2 

0.35 

0.75 

0.40 

0 

9 

Gibson 

IN 

38.4 

-87.8 

3.5 

56.5 

0.9 

17.5 

0.23 

0.90 

0.40 

2 

10 

Big Bend 

FL 

27.8 

-82.4 

4 . 2 * 

76.7 

2.3 

40.9 

0.53 

0.89 

0.31 

0 

11 

Laramie River 

WY 

42.1 

-104.9 

2.6 

41.8 

0.5 

19.7 

0.14 

0.66 

0.17 

2 

12 

Cumberland 

TN 

36.4 

-87.7 

4.4 

79.3 

1.3 

30.7 

0.23 

0.75 

0.32 

0 

13 

Roxboro 

NC 

36.5 

-79.1 

3.5 

67.2 

1.5 

29.8 

0.38 

0.92 

0.40 

0 

14 

Centralia 

WA 

46.8 

-122.9 

1.0* 

33.4 

0.6 

12.7 

0.23 

0.24 

0.40 

13 

15 

Navajo 

AZ 

36.9 

-111.4 

2.7 

28.5 

0.9 

24.1 

0.26 

0.88 

0.23 

0 

16 

Barry 

AL 

31.0 

-88.0 

2.0 

49.3 

0.5 

15.6 

0.26 

0.74 

0.17 

0 

17 

Powerton 

IL 

40.5 

-89.7 

2.1 

41.3 

0.9 

19.3 

0.21 

0.44 

0.36 

9 

18 

Mount Storm 

WV 

39.2 

-79.3 

3.3 

81.4 

1.6 

33.4 

0.41 

0.83 

0.51 

3 

19 

Dave Johnston 

WY 

42.8 

-105.8 

2.5 

62.3 

0.5 

24.3 

0.21 

0.76 

0.13 

8 

20 

Chalk Point 

MD 

38.5 

-76.7 

2.0* 

61.0 

2.6 

32.5 

1.08 

0.89 

0.66 

0 

21 

Wansley 

GA 

33.4 

-85.0 

2.5 

78.3 

1.2 

22.6 

0.37 

0.54 

0.33 

0 

22 

Allen Fossil 

TN 

35.1 

-90.1 

2.4 

82.1 

1.0 

20.0 

0.43 

0.88 

0.35 

0 

23 

Kingston 

TN 

35.9 

-84.5 

2.5 

87.6 

1.1 

24.9 

0.47 

0.86 

0.41 

11 

24 

JM Stuart 

OH 

38.6 

-83.7 

2.8 

66.0 

2.2 

35.4 

0.62 

0.73 

0.64 

6 

25 

Dolet Hills 

LA 

32.0 

-93.6 

1 . 8 * 

58.0 

0.1 

2.8 

0.14 

0.47 

0.17 

3 

26 

Johnsonville 

TN 

36.0 

-88.0 

1 . 8 * 

45.5 

1.2 

27.0 

0.51 

0.66 

0.31 

0 

27 

St. Johns River 

FL 

30.4 

-81.6 

2.6 

60.5 

1.9 

34.6 

0.53 

0.73 

0.31 

0 

28 

Clifty Creek 

IN 

38.7 

-85.4 

2.1 

53.4 

1.9 

30.5 

0.65 

0.57 

0.60 

4 

29 

Conemaugh + Homer City 

PA 

40.4 

-79.1 

2.9 

37.1 

2.1 

27.2 

0.49 

0.39 

0.78 

10 

30 

Hunter 

UT 

39.0 

- 111.0 

2.0 

38.2 

0.5 

14.6 

0.21 

0.60 

0.37 

7 

31 

Antelope Valley 

ND 

47.4 

-101.8 

1.6* 

44.0 

0.1 

4.8 

0.16 

0.32 

0.19 

15 

32 

Marshall 

NC 

35.6 

-81.0 

1 . 3 * 

35.5 

1.6 

26.5 

0.79 

0.59 

0.39 

0 

33 

Big Stone 

SD 

45.3 

-96.5 

u* 

39.2 

0.3 

11.5 

0.21 

0.73 

0.18 

23 

34 

Boswell 

MN 

47.3 

-93.7 

1 . 1 * 

50.4 

0.1 

6.4 

0.15 

0.45 

0.19 

22 

35 

George Neal North 

IA 

42.3 

-96.4 

00 

o 

31.3 

0.2 

7.7 

0.32 

0.39 

0.22 

11 

36 

Harllee Branch 

GA 

33.2 

-83.3 

1.7* 

35.2 

1.1 

25.0 

0.38 

0.61 

0.27 

0 

37 

Mill Creek 

KY 

38.1 

-85.9 

1 . 1 * 

40.5 

0.8 

12.9 

0.49 

0.30 

0.50 

0 

38 

Monticello 

TX 

33.1 

-95.0 

0 . 7 * 

20.0 

0.4 

11.2 

0.56 

0.44 

0.23 

0 

39 

Ghent 

KY 

38.7 

-85.0 

1 . 0 * 

39.6 

1.8 

27.4 

0.80 

0.31 

0.57 

2 

40 

Coronado 

AZ 

34.6 

-109.3 

0.3* 

10.1 

0.2 

9.3 

0.18 

0.34 

0.15 

0 

41 

Scherer 

GA 

33.1 

-83.8 

0 . 5 * 

10.9 

1.0 

22.0 

0.84 

0.43 

0.31 

0 

42 

Greene County 

AL 

32.6 

-87.8 

LD 

d 

28.4 

0.4 

13.8 

0.65 

0.65 

0.16 

0 

43 

Craig 

CO 

40.5 

-107.6 

0.6* 

22.3 

0.6 

25.0 

0.86 

0.89 

0.25 

30 

44 

Grand River Dam 

OK 

36.2 

-95.3 

0.3* 

8.4 

0.6 

14.3 

0.34 

0.26 

0.25 

0 

45 

Muskogee 

OK 

35.8 

-95.3 

0.4* 

9.4 

0.6 

17.3 

0.17 

0.10 

0.23 

0 

46 

White Bluff 

AR 

34.4 

-92.1 

-0.2* 

-3.9 

0.0 

0.5 

0.24 

0.22 

0.21 

1 

47 

Big Cajun2 

LA 

30.7 

-91.4 

0.1* 

2.3 

0.3 

6.7 

0.51 

0.16 

0.25 

0 

48 

Belle River + St. Clair 

MI 

42.8 

-82.5 

-0.1* 

-4.6 

1.1 

17.5 

-0.07 

0.00 

0.61 

21 

49 

Gadsden 

AL 

34.0 

-86.0 

0.2* 

42.1 

0.6 

16.5 

0.74 

0.16 

0.21 

0 

50 

LelandOlds + MiltonRYoung 

ND 

47.2 

-101.3 

-0.5* 

-10.4 

0.8 

21.9 

-0.21 

0.15 

0.24 

20 

51 

Rockport 

IN 

37.9 

-87.0 

-0.1* 

-2.9 

1.8 

27.7 

-0.91 

0.11 

0.42 

4 

52 

Kincaid 

IL 

39.6 

-89.5 

-0.6* 

-26.3 

0.8 

21.2 

0.00 

0.00 

0.35 

4 

53 

Four Corners + San Juan 

NM 

36.7 

-108.5 

0.0* 

-0.3 

1.0 

15.0 

-0.19 

0.18 

0.45 

2 

54 

Martin Lake 

TX 

32.3 

-94.6 

-0.7* 

-19.4 

0.3 

8.6 

-0.22 

0.06 

0.24 

0 

55 

Montour 

PA 

41.1 

-76.7 

0.6* 

25.6 

0.8 

14.4 

1.04 

0.94 

0.66 

10 


AE PP and ANOj are calculated as the mean of 2005 and 2006 minus the mean of 2010 and 2011 ; the total sum of emissions includes only those days where OMI data are 
available so that AE PP is less than the total change in a facility’s emissions over a given time period. 

p is the slope of the line fit to annual mean OMI data and annual total CEMS data in Fig. 2 and SI (right column); alternately, one could define p as the change in NO 2 and E pp 
between two specific years, such as 2005 and 201 1 , which gives a similar value for p as the slope of a linear fit to the data when AE PP is large. 
r 2 is the correlation of annual mean OMI data and annual total CEMS data (Fig. 2 and SI; right column). 

StdErr is the standard error of the mean of NO 2 . 

N < 3 is the number of months between 2005 and 201 1 which have less than three days to create the monthly average. 

We indicate those facilities where the regional background is large relative to the signal of the facility with an in the column containing AE PP . Our criteria for this 

determination is AE PP <2 kTon (see Fig. 6) and/or location within regionally polluted areas. 


for each site so that ANO 2 can be treated as linearly proportional to 
AE 1 (e.g., Martin et al., 2003; Kaynak et al„ 2009; Kim et al„ 2009). 

3. Results 

Over the Aura record, 2005—2011, NO* emissions from electric 
power generation decreased by 48% in the US (http://www.epa. 


gov/ttn/chief/trends/index.html). The CEMS data indicate that 
there was a large (>50%) decrease between 2005 and 2011 in 
annual emissions at 22 power plants that we include in our study 
(Table 1 ), presumably because of the implementation of new ECDs. 
Emissions at most of the facilities decreased by >20%, while 
emissions at seven facilities did not change or increased. At many 
facilities, emissions decreased rapidly, but they decreased 
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Fig. 2. (left) Monthly mean NOj (black line; x 10 15 molecules cm -2 ) and E pp (blue dotted line; kTon) data from 2005 to 2011 for four power plants. Vertical black lines represent the 
standard error of the mean of the OM1 data. The sample size (N) is the number of days with data used to create monthly means. The annual mean NOj data are represented with a 
green line and the monthly median data as an open red diamond, (right) Annual mean NO 2 ( x 10 15 molecules cur 2 ) versus annual E pp (kTon) data. The red numbers represent the 
years that correspond to the annual means (e.g., “09’' = 2009). The correlation (r 2 ) of the data is shown along with the slope (m) and y-intercept (b) of a line fit to the data, m is p as 
shown in Table 1. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


incrementally at others, which we know to be associated with, for 
instance, the implementation over time of ECDs on specific units 
within a facility. 

3.3. Response of OMI N0 2 to the implementation of ECDs 

Fig. 2 (left column) shows monthly total E pp and monthly 
mean NO 2 for several facilities over the OMI record (Figure SI 


shows this information for all facilities listed in Table 1.). The 
Crystal River facility (ID #1; Fig. 2a) in Florida had the largest 
AE PP (70%; Table 1). The CEMS data show that emissions began 
decreasing rapidly during the installation of ECDs that came 
online in June 2009 (Unit 5) and May 2010 (Unit 4). There were 
concomitant decreases in NOj with a 37% overall reduction 
(Table 1). The correlation of the annually mean NO 2 and annual 
total E pp (Fig. 2, right column) is high (r 2 = 0.91). p (i.e., the slope 
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All Power Plants 



Northern US 



Southern US 



A N0 2 t [x 10 ,5 molec/cm 2 ] 


^ [kTon] 

0.0 2.0 4.0 6.0 8.0 

Fig. 3. (top) AE PP (kTon) as compared to ANOj (x 10 15 molecules cirr 2 ) as the mean of 
2005 and 2006 minus the mean of 2010 and 2011 (Table 1). The colored dots indicate 
the magnitude of the mean E pp of 2005 and 2006. The number associated with each 
point corresponds to a particular power plant identified in Table 1. n is sample size (i.e., 
the number of power plants) used in the correlation statistic (r) and line fit, where m is 
the slope and b is the y-intercept. (middle) The same as (top), but for only those fa- 
cilities at latitudes >36.5° N. The horizontal lines represent the standard error of the 
means of the OMI data, (bottom) The same as (middle), but for only those facilities at 
latitudes <36.5° N. 


(m) of the line in Fig. 2 and SI; right column) at the Crystal River 
facility is 0.28. 

p at the Bowen facility (ID #8) in Georgia is similar (0.35) to the 
Crystal River facility, though the correlation is somewhat lower 
(r 2 = 0.75; Fig. 2b). ECDs were operated at Bowen during the ozone 
season through 2008, but year-round afterward (Fig. 2b; left col- 
umn). Overall, E pp decreased by 70% from 2005 to 2011 with a 


corresponding decrease in NOj of 30%. Although this facility is 
generally upwind of the Atlanta metropolitan area, the NO 2 column 
is likely influenced to some degree by this urban source, depending 
on meteorology and season, which may explain the scatter in E pp 
and NOj in Fig. 2b (right column). Nevertheless, the impact of year- 
round ECDs on NO 2 is clear from the beginning of 2009. 

We found that p’s are <1 at all but two of the facilities (Table 1 ). 
There is a wide range of values, but there is no clustering, such as 
with latitude. In the next section, we discuss sources of variation of 
p among the facilities. 

3.2. Sources of variation in the response 

Fig. 3 shows AE PP and ANO 2 (Table 1) for all power plants. 
Overall, the correlation (r 2 = 0.31) is weak and does not improve 
much when only facilities are considered where E pp > 4 kTon in 
2005. The poor correlation occurs whether the absolute or relative 
changes are considered. The correlation for the facilities in the 
southern US is better (r 2 = 0.58; n = 25) than in the northern US 
(r 2 = 0.18; n = 30; Fig. 3), though it is important to note that p’s at 
most individual facilities are generally scalar, including in the 
northern US (Figure SI ) (There are too few facilities in the western 
US with which to draw any conclusion about a possible systematic 
bias between eastern and western facilities.). We chose 36.5°N 
latitude (shown in Fig. 1 ) to separate the northern and southern US 
as the range of the standard error of the means of the OMI data for 
the individual facilities is 0.1— 0.4 below this latitude and 0.1— 0.8 
above this latitude (Fig. 3; Table 1 ). In the following subsections, we 
discuss the sources of variation of p for the power plants, including 
those that cause the differences between facilities in the northern 
and southern US. 

3.2.1. Magnitude of emissions reduction 

For most facilities, we found that p is scalar for all intents and 
purposes and ANO 2 is well correlated with AE PP given that AE PP was 
large (Fig. 1 ; Table 1 ), which is consistent with the findings of lu 
and Streets (2012). The correlations (r 2 ) between annual NO 2 and 
E pp (Table 1 ) are >0.5 at 32 of the 55 facilities and, not surprisingly, 
rise linearly with increasing AE PP , which will be discussed further in 
Section 3.2.6. 

As with all satellite data, it is important to consider the issue of 
the signal-to-noise ratio (SNR). For our purposes, this means that 
the SNR increases with the magnitude of E pp . Some of the facilities 
in Table 1 had relatively small annual emissions in 2005 so that 
meteorological variations and large changes in regional N0 2 levels, 
for instance, may obscure their p’s. Not all the power plants in our 
study used ECDs and some had relatively small variations in annual 
emissions. We included these facilities to help us understand what 
factors influence p. 

3.2.2. Retrieval issues 

Errors are introduced into the retrieval during the conversion of 
the measured OMI slant column to a more useful vertical column 
using a tropospheric AMF, a complex function of information on a 
priori NO 2 profile shapes, surface albedo, clouds, aerosols (not 
implicitly accounted for), etc. (e.g., Boersma et al., 2011). The use of 
coarsely-resolved retrieval parameters (e.g., N0 2 profile shapes, 
surface albedo) could introduce large errors in retrievals at places 
where these parameters have large spatial variability (Zhou et al., 
2009; Boersma et al., 2011), such as in mountainous and desert 
areas in the western US. For example, the emissions remained 
relatively stable over our study period at the Four Corners/San Juan 
facility (ID #53) in New Mexico, the facility with one of the highest 
annual emissions (Fig. 2c). Though new ECDs were not installed, 
the year-to-year variation in E pp was larger at this facility than AE PP 
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for some smaller facilities that implemented ECDs. However, the 
correlation (r 2 = 0.18) associated with p is weak, which may result 
from the parameters used in the AMF as the facility is located in the 
desert (i.e., high surface reflectivities) and near mountains char- 
acterized by variable snow cover. 

3.2.3. Statistical significance 

The number of individual days with OMI data (i.e., sample size 
(N) in Fig. 2 and SI) is typically < 10 month -1 , so that the standard 
error of the mean is oftentimes large. In these situations, the 
monthly average is not statistically significant. This issue is com- 
pounded for power plants at higher latitudes or elevations, such as 
the Big Stone (ID #33; Figure SI ) facility, as OMI data are filtered for 
snow cover. At this facility, 23 months have <3 days of data with 
which to create the monthly average (Table 1 ), so that the annual 
average is weighted more heavily to spring, summer, and fall than 
winter. It is worth noting that the correlation of p (r 2 = 0.73) is high 
for this facility. 

As N is not large, the monthly OMI data, as gridded for use in this 
study, may be skewed by outliers and bad data at all facilities. There 
are two winter months with obvious bad data, possibly because of 
improper filtering for snow and ice, at the Monroe facility (ID #7) in 
Michigan (Figure SI). However, the impact of this bad data is not 
obvious in the correlation (r 2 = 0.58) of p. Other facilities with 
suspect data during winter include, for instance, Gibson (ID #9), 
Boswell (ID #34), and Four Corners/San Juan (ID #53). 

3.2.4. Proximity to urban sources 

At the Big Bend facility (ID #10) near Tampa, Florida, ECDs were 
brought online in 2008 (Unit 3), 2009 (Unit 2), and 2010 (Unit 1), 
decreasing emissions by 77%. Similar to p (0.28) at the Crystal River 
facility (ID #1), which is also in Florida, p (0.53) at the Big Bend 
facility is scalar, but twice as high; NO 2 and E pp are well correlated 
(r 2 = 0.89). ANO 2 at the Crystal River facility is smaller than at the 
Big Bend facility despite AE PP being larger for the Crystal River fa- 
cility. Due to proximity, the urban plume of Tampa influenced NO 2 
at the Big Bend facility (not shown), particularly in the earlier years 
of our study period. From 2005 to 2011, the OMI data indicate that 
NO 2 over Tampa decreased by more than 50% (~2.5 x 10 15 mole- 
cules cm -2 ), which, coupled with the large AE PP , explains the larger 
p as compared to the Crystal River facility. That is, p for the Big Bend 
facility is convolved with the large decrease of NO 2 in the urban 
plume of Tampa (i.e., ANO-?). For facilities near large emitters, 
including cities, the OMI data could be filtered by wind direction to 
minimize the influence of these other sources. 

3.2.5. Seasonal variation 

The influence of the seasonal cycle in NO 2 associated with var- 
iations in temperature and sunlight is readily apparent in Fig. 2b at 
the Bowen facility (ID #8) and at numerous other facilities 
(Figure SI). At facilities, such as Paradise (ID #3), New Madrid (ID 
#4), and Gibson (ID #9), the seasonal cycles in NOj continue even 
after the ECDs were routinely used year-round. It is worth noting 
that the correlation of monthly NO 2 and E pp may be artificially 
enhanced by the coincidence of the seasonal minimum of the 
chemical lifetime of N0 2 in summer and the use of ECDs during the 
ozone season only (e.g., the Montour facility (ID #55); Fig. 2d). 

As discussed in Section 2.1, we use all available data, regardless 
of season, to calculate p as many facilities operated ECDs in sum- 
mer during our entire study period. The advantage of using all data 
is that the sample size (JV) is larger, thus improving statistical 
significance as discussed in Section 3.2.3. Using all available OMI 
data will cause variation in the p’s calculated for the facilities. The 
seasonal variation of the chemical lifetime is greatest at higher 
latitudes, which partly explains why ANO 2 and AE PP from the 
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Fig. 4. Percent change of the annual mean OMI NOj (dashed lines) and N0 2 from AQS 
surface sites (solid lines) relative to 2005 for the whole US ("All US") and four quad- 
rants (“NE" = northeast”; “SE” = southeast; “NW” = northwest; “SW” = southwest). In 
total, 517 AQS sites are included. The OMI data were sampled for the 0.1° x 0.1° 
gridboxes in which the AQS sites lay. We used all hourly AQS data to estimate the 
percent change. 


individual facilities are better correlated in the southern (r 2 = 0.58) 
than in the northern US (r 2 = 0.18; Fig. 3). However, the calculation 
of p for an individual facility should not be adversely impacted as a 
similar distribution of data over the course of a year is used for all 
years. 

To understand the seasonal variability, we calculated p for each 
of the four seasons for each facility. In general, there is significant 
variability (>50%) in the seasonal p’s for the typical facility, 
particularly ones at higher latitudes. At the Crystal River facility (ID 
#1 ) in Florida, the seasonal p’s are similar (i.e., within ~ 30%), which 
is not surprising given the plant’s southerly location. The seasonal 
p’s show more variation (~50%) at the nearby Big Bend facility (ID 
#10), though this facility is impacted by the urban plume of Tampa 
as discussed in Section 3.2.4. On the other hand, there is consid- 
erably more variability in the seasonal p’s at the Cardinal/W. H. 
Sammis facility (ID #2) in Ohio and the New Madrid facility (ID #4) 
in Missouri, which are both located at higher latitudes than the 
Crystal River facility and in areas with higher regional NO 2 levels. 

In general, the seasonal p’s for spring, summer and fall tend to be 
more similar for a typical facility with the seasonal p for winter 
being the outlier. An exception is that the seasonal p’s for summer 
are less meaningful at facilities in which ECDs were used during the 
ozone season over the entire study period because of the low 
emissions and, subsequently, low SNR of the OMI data. In addition, 
the magnitude of NO 2 is seasonally lowest in summer as the 
chemical lifetime is seasonally shortest. One would expect that 
seasonal p’s for spring and fall are similar at a given facility because 
of the similar chemical lifetimes in these two seasons. Generally, 
this is the case, particularly for the high-emitting facilities. 

The seasonal p’s for winter tend to show considerable variability 
because of the latitudinal-dependence of the seasonal variation of 
the chemical lifetime and because of missing data as discussed in 
Section 3.2.3. In addition, the regional NO 2 levels are seasonally 
highest in winter because the chemical lifetime is seasonally 
longest, so that the ratio of NO 2 from the power plant relative to the 
regional level is seasonally lowest at many facilities. 

The methodology that we present in this manuscript to assess 
the impact of the implementation of an ECD on the N0 2 level above 
a facility may be tailored for specific applications. For instance, one 
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Fig. 5. Annual mean NOj (x 10' 5 molecules cnr 2 ) for a) 2005 and b) 2011 for two re- 
gions: the Ohio River Valley (OHRV), which has a high concentration of power plants, and 
the densely populated Northeast Corridor (DC-NYC). Data <2 x 10 15 molecules cm~ 2 , 
which have low SNRs, are shown as white in these two regions. The c) absolute and d) 
percent changes in the mean data for the two regions relative to levels in 2005. 


may decide to exclude the winter season from the analysis because 
of issues associated with NO* lifetime and statistical significance as 
discussed above. We repeated our analysis, excluding the winter 
season, and found that the conclusions of our study remain un- 
changed. The correlations of AE PP and AN0 2 changed modestly for 
most facilities, though we did not find an improvement in the 
scatter observed in Fig. 3 as indicated by the correlation statistic (r). 
r 2 for the northern US decreased from 0.18 to 0.07 and remained the 
same for the southern US. Based on the discussion in this section, it 
is important to understand that the value of p will depend on the 
choice of months used in the analysis, particularly in regions with 
seasonal variation in the NO* lifetime. 


Fig. 6. p (xlO 15 molec (cm 2 kTon) 1 ) versus AE PP (kTon) for the individual facilities. 
The vertical dashed line separates AE PP into two categories: AE PP > 2 kTon, 
AE PP < 2 kTon. 

3.2.6. Variations in regional JV0 2 levels (ANO 2 ) 

Annual mean surface concentrations of N0 2 decreased 33% na- 
tionally between 2001 and 2010 (EPA, 2012), primarily from de- 
creases in emissions from fuel combustion in electrical utilities and 
vehicles (e.g., McDonald et al., 2012). From 2005 to 2011, the com- 
bined NO* emissions reduction from electrical utilities and vehicles 
was about 37% (http://www.epa.gov/ttn/chief/trends/index.html) 
with two-thirds of the decrease being attributed to the reduction in 
the mobile source. The OMI data confirm that regional N0 2 levels 
decreased substantially in many areas of the US during our study 
period (e.g., Russell et al., 2012). Fig. 4 shows the percent change 
relative to 2005 of EPA’s Air Quality Monitoring System (AQS; http:// 
www.epa.gov/ttn/airs/airsaqs/) N0 2 data and the corresponding 
OMI data above the AQS stations. The data are averaged over the 
whole US and over four quadrants. The reductions by 2009 range 
from 20 to 30% for the AQS sites as grouped in the quadrants, but 35— 
40% for OMI data; the northwest quadrant is an outlier. However, the 
overall shapes of the trends in both datasets are similar. The 
discrepancy in the magnitudes of the trends of the AQS and OMI data 
may occur as the OMI detects changes in N0 2 throughout the whole 
troposphere, while monitors at the AQS sites sample near-surface 
air. Thus, the OMI detects the reductions in N0 2 associated with 
both mobile and power plant sources, while the AQS surface mon- 
itors preferentially sample reductions in mobile sources as power 
plant plumes are located aloft predominately. 

Though the chemical lifetime of NO* is relatively short, NO„ 
emissions upwind influence N0 2 columns downwind (e.g., Turner 
et al., 2012), which can lead to elevated regional N0 2 levels (As 
discussed in Section 2.2 for our analysis, we selected facilities least 
affected by other sources within a 0.4° longitude x 0.4° latitude 
area around the facility (i.e.,/ p0 wer > 0.90)). Over our study period, 
some of the largest changes in regional levels occurred in the 
heavily populated region extending from Washington, DC to New 
York City (i.e., the Northeast Corridor), and the industrialized Ohio 
River Valley, where eight of the power plants selected for this study 
are located; the Chalk Point facility (ID #20) is the only facility that 
met our criterion (i.e.,/ p0 wer > 0.90) for selection in the Northeast 
Corridor. Fig. 5 shows that regional N0 2 levels decreased by 30— 
40% from 2005 to 2011 in both of these regions, though the absolute 
decrease was much higher in the Northeast Corridor. Most of the 
power plants in Table 1 are located in areas with lower regional N0 2 
levels in 2005 than in the Northeast Corridor and Ohio River Valley. 
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Fig. 6 shows the relationship between AE PP and p. For power 
plants with AE PP > 2 kTon, AE PP is generally large relative to the 
change in the regional NO 2 level, p’s for these facilities are between 
0.12 and 0.62 with a mean of 0.36. For facilities with AE PP < 2 kTon, 
there is a wider range of p’s (i.e., between -0.91 and 1.08), indicating 
that a change in the regional NO 2 level, if large, can influence p in a 
non-negligible way. We attempted to find a method for removing 
the influence of a change in the regional NO 2 level in a general way 
applicable to all facilities. However, we found that the change in the 
regional NO 2 level can vary widely (e.g., with meteorological vari- 
ability), requiring careful processing of the data for each facility. As 
an example, NOj’s were high at the nearby facilities of White Bluff 
(ID #46) and Dolet Hills (ID # 25) in the winters of 2009—10 and 
2010—11 (Figure SI), which we found to be caused by stagnant 
meteorological conditions that allowed regional NO 2 levels to build. 

4. Summary 

We conclude that it is practical to use OMI N0 2 tropospheric 
column data to assess changes of emissions from power plants that 
are associated with the implementation of emission control devices 
(ECDs), though careful interpretation of the data is necessary. We 
showed that there is a clear response for OMI NO 2 data to NO* 
emission reductions from power plants associated with the 
implementation of ECDs on both monthly and annual timescales. 
This response is scalar for all intents and purposes, whether the 
reduction is rapid or incremental over several years. However, the 
response is variable among the power plants, even those with the 
greatest absolute decrease in emissions. We discussed some of the 
causes of this variability, which include the magnitude of a facility’s 
NO* emissions, seasonal variation of the NO* lifetime, proximity to 
urban areas, changes in the regional NO 2 levels, lack of statistical 
significance, and retrieval issues. Ideally, one should use a CTM to 
account for several of these causes of variability, though this would 
limit the practical application of space-based data for air quality 
purposes because of computational expense. However, we show 
that this step is not necessary if the change in the facility’s NO* 
emissions is large. 

Using space-based N0 2 columns to assess changes in power plant 
NO* emissions will likely become more quantitative as the OMI 
retrieval procedure continues to evolve, such as through the use of 
improved and finely-resolved information of surface parameters. In 
addition, two planned sensors promise enhanced capabilities as 
compared to OMI: i) the European Space Agency Tropospheric 
Ozone Monitoring Instrument (TROPOMI; http://www.knmi.nl/ 
samenw/tropomi/lnstrument/), an OMI follow-on instrument with 
finer horizontal resolution, and ii) the NASA Tropospheric Emis- 
sions: Monitoring of Pollution (TEMPO; http://science.nasa.gov/ 
missions/tempo/) instrument, an OMI-like instrument that will be 
in geostationary orbit, collecting data throughout the day as 
opposed to one overpass per day as with OMI. 

Many of the facilities included in our study were already using 
ECDs during the ozone season before the start of the OMI data re- 
cord. A next step would be to repeat our analysis over the SCIA- 
MACHY, GOME, and GOME-2 records, similar to the study of Lu and 
Streets (2012), extending the period of study to 1996. The limitation 
of this approach is that the horizontal resolutions of data from 
these instruments are coarser than OMI data, which would make it 
more difficult to isolate the signal of individual facilities from sig- 
nals of other nearby sources. 
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